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Interference alignment is a signaling technique that provides high multiplexing gain in the interference 
f — . channel. It can be extended to multi-hop interference channels, where relays aid transmission between 

I— I sources and destinations. In addition to coverage extension and capacity enhancement, relays increase the 

H 

l-H multiplexing gain in the interference channel. In this paper, three cooperative algorithms are proposed 

O for a multiple-antenna ampUfy-and-forward (AF) relay interference channel. The algorithms design the 

transmitters and relays so that interference at the receivers can be aligned and canceled. The first algorithm 
^ minimizes the sum power of enhanced noise from the relays and interference at the receivers. The second 

and third algorithms rely on a connection between mean square error and mutual information to solve 

in 

the end-to-end sum-rate maximization problem with either equality or inequality power constraints via 
matrix-weighted sum mean square error minimization. Since we can find a globally optimal solution 
^ I in each iteration, the resulting iterative algorithms are convergent. Simulations show that the proposed 

^ algorithms achieve higher end-to-end sum-rates and multiplexing gains that existing strategies for AF 

relays, decode-and-forward relays, and direct transmission. The first algorithm outperforms the other 

H 

algorithms at high signal-to-noise ratio (SNR) but performs worse than them at low SNR. Thanks to 
power control, the third algorithm outperforms the second algorithm at the cost of additional overhead. 
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I. Introduction 

Relay interference channels model networks where a stage of intermediate nodes, called relays, help 
multiple transmitters communicate with their receivers using shared radio resources |[T|-||4|. Upcoming 
cellular standards are considering relay communication for coverage extension and capacity enhance- 
ment Q, Q. Prior work, however, shows that single- antenna relays do not work well in the presence 
of co-channel interference ||7|, ||8|. In this paper, we consider multiple-antenna relay systems to take 
the advantage of the interference management capability of multiple-input multiple-output (MIMO) 
communication. Many interference management strategies have been proposed for the MIMO single- 
hop interference channel ||9|-|[T3|. Although these single-hop results can be applied separately for the 
transmitter-relay hop and for the relay-receiver hop, even higher sum-rates can be achieved if the relays 
are configured jointly |14|, p3| . Obtaining the most from relay interference channels requires advanced 
interference management strategies that jointly configure the transmitters, relays, and receivers. 

A general challenge to designing algorithms for the interference channel is that the sum capacity is 
unknown. The multiplexing gain of a network is a first-order approximation of its sum-capacity at high 



signal-to-noise ratio (SNR) |16|. Interference alignment is a multiplexing gain maximizing signaling 



technique for the single-hop interference channel |13|, achieving the maximum number of degrees of 
freedom. The idea is to arrange the transmitted signals such that interference is constrained within only a 
portion of the signal space observed by each receiver, leaving the remaining portion for interference-free 
detection of the desired signal |9|. The maximum multiplexing gains achievable through interference 
alignment, however, depend on the characteristics of the interference channels. For a symmetric MIMO 
interference channel with constant channel coefficients, the maximum multiplexing gain is upper-bounded 
by the total number of antennas at a transmitter-receiver pair regardless the number of pairs |11|, |[T2j. 
Note that the bound is tight in certain cases and corresponds to the total available spatial dimensions of 
a pair. Increasing the number of spatial dimensions in the network, using for example relays, is one way 
to improve the maximum achievable multiplexing gain. 

Relays can be classified based on their signal processing operation, among which the most popular are 
decode-and-forward (DF - the relays decode the received signals then re-encode before retransmitting) 
and amplify-and-forward (AF - the relays apply linear signal processing to the signal before forwarding). 
Without decoding the received signals, AF relays need no knowledge of the codebooks used by the 
transmitters and Ukely have lower baseband complexity and fast signal processing. In addition, transparent 
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to the modulation and coding of the signals, AF relays are more suitable for applications in heterogeneous 
networks comprising many nodes of different complexity or even standards | [T7[ . In this paper, we focus 
on a half-duplex MIMO AF relay interference channel. Since half-duplex relays cannot transmit and 
receive at the same time, they are more practical than full-duplex relays. 

Several interference management strategies designed specifically for the one-way AF relay interference 



channel have been proposed |15|, |18|-|32|. Although relays cannot improve the multiplexing gains of 
the single-antenna fully-connected interference channel with time- varying or frequency-selective channel 
coefficients | |T8[ , they are beneficial for reducing the number of independent channel extensions needed 
to align interference at the receivers [19|. Prior work often considers networks operating in special 



circumstances. It is assumed in |15|, |20|-|24| that there are enough antennas at the relays to cancel 
all interference on the reception and then to nullify all interference on the retransmission, allowing 
multiplexing gains to scale linearly with the number of users. Other prior work considers only small 



networks with up to three pairs to derive some kind of closed-form strategies |25 |-p7|. Prior work 



in 1 28 1 considers design problems with different objective functions including sum power minimization 



and minimum SINR maximization. In addition, the algorithms in |28| are applicable only for single- 



antenna receivers. Prior work in |30|-|32| develops noncooperative resource allocation strategies for AF 
relay networks while our work focuses on centralized algorithms for cooperative resource allocation. In 
this paper, we consider a general setting in the sense that we assume no special constraints on the number 
of wireless nodes or the number of antennas at a node. The closest AF relay model to ours is considered 



in 1 29 1, which is only for single- antenna transmitters and receivers. Further, they assume no crosslinks 
from relays to receivers, resulting in an oversimplified design problem. 

Sum-rate maximization problems are nonconvex and NP-hard, i.e., their global optima cannot be found 
in a polynomial time. It is even challenging to find their good local optima corresponding to interference 
aligned solutions using gradient-based algorithms from arbitrary initializations because those solutions 



have very narrow regions of attraction |33 1. Thus, we propose to formulate three new design problems that 
have exactly the same constraints with sum-rate maximization problems but with better-behaved objective 
functions to find high-quality solutions in terms of sum-rate maximization. Based on the observation that 



interference alignment solutions make total leakage power go to zero |19|, |34|, |35|, we formulate 
a problem that aims at minimizing the sum power of interference and enhanced noise from the relays. 
Based on a relationship between achievable rates and mean squared error (MSE), we formulate two matrix- 
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weighted sum-MSE minimization problems, either without power control or with power control. The key 
is that they have the same stationary points as their corresponding end-to-end sum-rate maximization 
problems. Although the newly formulated optimization problems are still nonconvex and NP-hard, they 



may be easier to solve. Next, we propose to adopt an alternating minimization approach |36|, |37| 
to develop iterative algorithms for solving the newly formulated problems. In each iteration, all but 
one variable is fixed and we focus on designing the remaining variable by solving a single- variable 
optimization problem obtained from the corresponding original multi-variable problems. Since we are 
able to find a global optimum for the single-variable optimization problem in each iteration, the proposed 
algorithms are guaranteed to converge. Note that the power constraints at the relays depend on both the 
transmit precoders at the transmitters and the processing matrices at the relays, adding more constraints to 
the design problems. Thus, it is not straightforward to extend the methods used for the single-hop design 
problems to solve the two-hop design problems. Our initial results in this paper were reported in |38|. 



Compared with |38|, this paper presents three different algorithms, has more discussion of convergence 
and provides simulations that emphasize the achievable end-to-end sum-rates and multiplexing gains. 

We use Monte Carlo simulation to evaluate the average end-to-end sum-rates and multiplexing gains 
achievable through the proposed algorithms. First, the numerical results confirm the convergence of the 
proposed algorithms as expected. Second, over the iterations of the total leakage minimization algorithm, 
the true interference dominates at the beginning but is canceled quickly; after that, the enhanced noise 
from the relays becomes dominant. This means that relay-aided interference alignment should take into 
account the enhanced noise from the relays. Third, the total leakage minimization algorithm achieves lower 
end-to-end sum-rates than the others at low-to-medium SNR values because it ignores the desired signal 
power and noise power at the receivers. Nevertheless, the MSE-based algorithms result in unfairness, 
i.e., some users have much smaller rates than the others. Thus, the MSE-based algorithms achieve lower 
end-to-end sum-rates and multiplexing gains than the total leakage minimization algorithm at high SNR. 
One reason for this is that the MSE-based algorithms may either turn off some data streams or nullify 
the desired signals to some receivers. Fourth, for fixed numbers of antennas at the transmitters and at the 
receivers, even with half-duplex loss, AF relays can provide larger end-to-end multiplexing gains than DF 
relays or direct transmissions. Finally, the results show that AF relays provide larger average achievable 
end-to-end sum-rates than do DF relays. The proposed algorithms also provide higher achievable end- 
to-end sum-rates than the existing AF relaying strategies that do not align interference at the receivers. 
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The organization of the remainder of this paper is as follows. Section [IT] describes the system model. 
Section llll] formulates the end-to-end sum-rate maximization problems and presents our proposed ap- 



proach. Section IV develops three cooperative algorithms that aim at finding high-quality solutions of the 
sum-rate maximization problems. Section |V] evaluates numerically the proposed algorithms. Section VI 
concludes this paper and suggests future research. 

Notation: We use normal letters (e.g., a) for scalars, lowercase and uppercase boldface letters (e.g., h 
and H) for column vectors and matrices. Ijv and 0^ are the identity matrix and all-zero matrices of size 
N X N. t'^j^(A) gives the eigenvectors corresponding to the n smallest eigenvalues of A. For a matrix 
A, A^ is the transpose matrix, ||A|||, the Frobenious norm. A* the conjugate transpose, and tr(A) 
the trace. vec(A) denotes the vec operator to transform A into a while vec~^(a) denotes the inverse 
operator, (g) is the Kronecker product. E[-] is the statistical expectation operator. denotes iteration 
index. ()t is used for transmitters' parameters, ()r for receivers', and ()x for relays'. 



II. System Model 

Consider a relay interference channel where M half-duplex AF relays aid the one-way communication 
between K pairs of transmitters and receivers, as illustrated in Fig. [T] Each transmitter has data for only 
one receiver and each receiver is served by only one transmitter. Each pair is assigned a unique index 
/c G /C = {1, • • • , K}. Transmitter k has N^^f, antennas while receiver k has N^^j. antennas for k ^ JC. 
Similarly, each relay is assigned a unique index m G 7W = {1, • • • , M}. Relay m has A'x,™ antennas 
for m G A^. The half-duplex relays cannot transmit and receive at the same time, thus the transmission 
procedure consists of two stages. In the first stage, the transmitters send data to the relays. In the second 
stage, the relays apply linear processing to the received signals and forward to the receivers. We assume 
the direct channels between the transmitters and the receivers are ignored by the second-stage receivers. 




Fig. L A relay interference cliannel wiiere M iialf-duplex AF relays aid the communication of K transmitter-receiver pairs. 
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We denote ilm,k £ £,^x,^xNT.k the matrix channel from transmitter k to relay m and Gk,m S 
(j~<AfR,fcX7Vx,™ lY^Q matrix channel from relay m to receiver k for k £ K. and m G ^A. We assume that 
perfect and instantaneous knowledge of iim,k and Gk,-m for e /C and m G is available at a central 
processing unit. Although this is a strict requirement, our results are still valuable since they show the 
substantial gains that can be achieved through coordination. Our results can be used as a benchmark for 
future work that makes more practical CSI assumptions. 

Let Sfc G ([^d,kXi jj^g transmit symbol vector at transmitter k, where < minjA'T.fc, -^R,fc} is the 
number of data streams from transmitter k to receiver A; for A; G /C. The transmit symbols are independent 
identically distributed (i.i.d.) such that E(s/;S^) = 1^^. Transmitter k uses a linear transmit precoder 
Ffc G c^T.fcXdfc jj^^p transmit antennas. Let p'^f^ be the maximum transmit power. The actual 

transmit power at transmitter k is pT,fc = tr(F|.Ffc). Let nx,m be spatially white, additive Gaussian noise 
at relay m with covariance E(nx,mnx m) — nJ-Nx ™ for m G M. With perfect synchronization, relay 
m observes the following signal 

K 

yx,m = ^ Hm,fcFfc Sfc + nx,m- (1) 

Let Urn G be the processing matrix at relay m. The transmit signal at relay m is given by 

K 

XX,m = UmYX,™ = ^ '^m%m,k^k + Umnx,m- (2) 

fc=l 

Relay m actually uses the following transmit power 

K 

)+^inMVmU*J. (3) 

There are two possible types of power constraints at the relays: i) a set of individual power constraints 
at the relays and ii) a sum power constraint at all the relays. Individual relay power constraints are often 



considered in the cellular system literature |39|, ||40|. While a sum power constraint is often considered 



in the ad hoc network literature to extend the lifetime of battery-powered relays |41 1, |42|. Let Pxm 
the maximum transmit power at relay m and p^^^ be the maximum sum transmit power at all the relays. 
When power control is considered, the individual relay power constraints are px,m < Pxm'^"^ ^ 
whereas the sum relay power constraint is X]m=i Px,m < p^^^. Without power control, the inequalities 
in the power constraint expressions are replaced by equalities. The following sections focus on the 



sum power constraint at the relays. Section |IV-D| discusses the applicability of individual relay power 
constraints while Section V simulates the impact of individual relay power constraints on achievable 
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end-to-end sum-rates. 

Let be spatially white, additive Gaussian noise at receiver k with covariance E(nFi,,fcn^ ^) 
'^Rfc-'-^Rfc- denote Qk,m = Gk^nJJm- Receiver k observes the following signal 



M 



Yk = ^ Gfc,m.Xx,m + nR^fc 



m=l 
K M 



M 



nx,m + nR,A 

g=l m=l 



(4) 
(5) 



771=1 



where 7"fc,g is the effective end-to-end channel from transmitter q to receiver k for k,q £ }C. Applying 
a linear receive filter e c^R.^xfifc to y^, receiver k obtains 



A/ 



yk = W*kT k,kSk + ^IT k,qSq + ^I0k,mnx,m + W^HR^fc . 



(6) 



desired signal q=^k 



777=1 



local noise 



enhanced noise from relays 



interference 

The pre-processing interference-plus-noise covariance matrix at receiver k is 

K M 

= Y '^k,qTl^g + Y '^X,m^k,mGk,m + "^R.fcld*.- 

9=1 



(7) 



7n=l 



q=^k 



For notational convenience, we denote {F} = {Ffcjf^p {U} = {lJm}^=i and {W} = {W^jf^^. 
We also denote Vlm,k — ^miim,k and yVk,m — W^G^ m for /c G /C and m G 7W. Table |l] summarizes 
the notation of equivalent channel gain matrices used in the paper for A;, g G /C and m £ A4. 

TABLE I 

Notation of equivalent channel matrices used in the paper for k,q £ K. and m e M. 



Equivalent channel matrix 


Definition 


'^7n,fc 


Hm,fcFfc 


Gk,m 




^m,k 








Tk,q 


5^rrt=l Q k,rri^'tn,q — 5^rn=l Gfc^mUmHjnfjrFg 



III. Problem Formulation and Proposed Approach 



We formulate the end-to-end sum-rate maximization problem in Section III-A and propose an approach 



to solving it in Section III-B 
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A. End-to-end Sum-Rate Maximization 

For tractable analysis, we assume Gaussian signaling is used. For a given {F} and {U}, the achievable 
rate for the k-th transmitter-receiver pair is maximized by using the linear MMSE receive filter | |43| 

^MMSE ^ ^rk,krik + ^k)-'rk,k, (S) 

where is given in Thus, we only need to focus on the design of {F} and {U}. Note that the 
corresponding maximum achievable rate is given by [43 1 

i?fc({F},{U}) = Iog2det(lrf,+n,,R^in,fc). (9) 

The sum of the end-to-end achievable rates is defined as 

K 

i?3um({F},{U}) = -j;iog2det(ErsE({p|^|U})). (10) 

k=l 

The end-to-end sum-rate maximization problem without power control is formulated as follows 

(OP-noPC): min -i?3um({F}, {U}) 
{F},{U} 

s.t. PT,fc=p¥,r,fc = i,--- (11) 

M 

J2P^^rn=pT''- (12) 

m=l 

Remark 1: The counterpart problem with power control can be obtained by replacing the equalities in 
the constraints by the inequalities. We denote it as (CP-PC). Power control may improve the end-to-end 
sum-rates at the expense of additional overhead because the central unit needs to inform the transmitters 
about both the norm and the shape of designed transmit precoders. 

Remark 2: The following (Fo,fe, Wo^/t, Uo,m) for A; E /C and m G satisfies the constraints of both 
(OP-PC) and (OP-noPC) 



Fo,fc = J^lN^,,xd,,keJC, (13) 
Wo,fc = y^l7VK..xd.,A: G/C, (14) 

Uo,m = y/apf^lN^^xN^,r^,rn £ M, (15) 

where a = ( Ef=i =i tr{Id,.N.^,n*^,k^^,klN.^.xdJ + Z^=i ^^x,m4,n^; 

Remark 3: (OP-PC) and (OP-noPC) are nonconvex and NP-hard. Moreover, even the smallest 
configuration of the MIMO AF relay interference channel with K = M = 2 and A^t = -Ax = A^r = 2 
requires the determination of twelve complex variables for the transmit precoders and relay processing 
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matrices, which makes even a brute force approach challenging. 
B. Proposed Approach 

Instead of directly solving for the globally optimal solutions of (CP-PC) and (OP-noPC), we aim 
at finding their high-quality solutions with reasonable computational complexity. To do this, in Section 



III-B 1 and Section III-B2[ we formulate two classes of new optimizations problems that have exactly the 



same constraints as (OP-PC) or (CP-noPC) but with different objective functions. 

1 ) Total Leakage Minimization: This section presents an approach for interference alignment in the 
AF relay interference channel, which is inspired by those for the single-hop interference channel in fl9\, 



1 34 1, |35|. The underlying observation for this approach is that when interference alignment is feasible, 
the sum power of the interference at all the receivers, also known as the leakage, is zero. From Q, 
there are three groups of unwanted signals at each receiver: i) interference, ii) enhanced noise from the 
relays, and iii) local noise. We denote X({F}, {U}, {W}) as the total leakage power of the AF relay 
interference channel. By evaluating the expectation and exploiting the independence of transmit signals 
Sfc for G /C and using the equality ||A|||. = tr(AA*), we obtain 

K K 

X({F},{U},{W}) = 5]5]tr(Wm,,n,gWfc). (16) 

fc=l q=l 
q^k 

In our opinion, the high SNR regime of the relay interference channel corresponds to high transmit 
power at both the transmitters and the relays. As a result, in addition to eliminating completely interfer- 
ence, we also need to eliminate the enhanced relay noise; otherwise, the enhanced relay noise power scales 
with the desired signal power, preventing the system from achieving high multiplexing gain. We denote 
J\f{{\J}, {W}) as the sum power of enhanced noise from the relays. By evaluating the expectations and 
exploiting the independence of the noise vectors at the relays, we obtain 

K M 

Ar({U},{W}) = Y.T.''lm^'("^kOk,mQlmWk). (17) 

k=l m=l 

Note that scaling down transmit power at either the transmitters or relays decreases the total leakage 
power at the receivers. For example, if ({(l/a)F}, {U}, {W}) is used instead of ({F}, {U}, {W}) 
where {(l/a)F} = {(l/a)Fi,--- ,(l/a)Fi^} and a > 1, then both the actual transmit power at the 
transmitters and the total leakage power decrease > 1 times. Thus, equality power constraints at 
the transmitters and relays are required to obtain a meaningful design problem. This means that power 
control should not be considered in the context of total leakage power minimization. In other words, 
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we do not use the total leakage minimization approach to find solutions to (CP-PC). Also to obtain a 
meaningful design problem, we add the orthonormal constraints on as W^W^ = I^^, for k ^ JC. 
Without such constraints, we can always use zero matrices as the receive filters to get zero total leakage 
power. Consequently, to find high-quality solutions of (CP-noPC), we propose to solve the following 
problem 

{TC): ^p^min^j X({F},{U},{W})+AA({U},{W}) 

s.t. iJT,fe=p¥,r,A:E/C, (18) 

M 

^PX,n.=pT'', (19) 
m=l 

WlWk = Id,,kGlC. (20) 

Note that (TC) is nonconvex and in general is NP-hard. Also, (TC) does not take into account the 
desired signal power and local noise at the receivers. 

Remark 4: The total leakage minimization problem formulated in (29] for an AF relay network is a 



simplified version of (TC). It is assumed in |29| that the transmitters and receivers are equipped with 
a single antenna. Each pair is aided by a dedicated multiple-antenna AF relay. The formulation in [29J 
does not consider power constraints at the relays. In addition, it is assumed that there are no cross-Unks 
for the transmissions from relays to receivers, i.e. Gk,q = for all k,q G IC and /c / (7. As a result, for 



fixed {F} and {W}, the algorithm in |29| can determine each Um separately. 



2 ) Sum Mean Squared Error Minimization: The section presents another approach that is based on a 



relationship between the achievable rates and MSE values at the receivers with Gaussian signaling |44|. 
This is inspired by prior work on the MIMO broadcast channel |j45j|, MIMO interference channel f33l, 
MIMO interference broadcast channel |(46l, and two-way relay channel ||47l, ||48|. Let Efc({E}, {U}, W^) 



be the MSE matrix at receiver k. After some manipulation, we obtain 

E,({E}, {u}, w,) = w^(n,fcn,, + Rfc)Wfc - wirk,k - n,fcWfc + 1^,. (21) 

Note that the MSE at receiver k, defined as MSE^ = tr(Efc({E}, {U}, W^)), is minimized by the 
linear MMSE receive filter W^'^^^. Moreover, it is well-estabUshed that |33|, |44|-|49| 

i?fc({F},{U}) = -log2det(Efc({F},{U},W^^^^SE))_ ^^l) 

We introduce auxiliary weight matrix variables {V} = (Vi, • • • , V^) that are square (V^. G C'^iX'^fc) 
and positive semidefinite for k £ K,. The weight matrices are just auxiliary variables for the 
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optimization technique and have no actual physical meaning. Define the matrix-weighted sum of MSE 
values as follows 

VFM5i?sum({F}, {U}, {W}, {V}) = J] ( tr (VfeEfc({F}, {U}, {W})) - logs det (V^) J . (23) 

k=i ^ ^ 
Then, we formulate the following weighted sum-MSE minimization problem 

iWMSe-noPC) : mill VFMS^sum({F}, {U}, {W}, {V}) 

{F},{U},{W},{V} VL J L J L J L J/ 

s.t. PT,fc=p¥,r'^e^ (24) 

M 

^PX,m=pT''- (25) 
m=l 

Similarly, we formulate (WM-SS-PC) by replacing the equalities in {WAdSS-noFC). Using the same 



steps in |46|, |49|, we can show that (WMSS-noPC) and (CP-noPC) have exactly the same stationary 
points if we use the linear MMSE receivers and choose the following matrix weights 

V°P*({F},{U},W^^^^SE) = E,i({F},{U},W^^^SE)^ (26) 

= ^d, + Tl^k^'f^'^Tk,k- (27) 

This observation is also true for [WMSS-PC) and (OV-PC). Thus, instead of directly solving (CP-noPC) 
(or (CP-PC)), we can focus on finding high-quality solutions to its corresponding weighted sum-MSE 
minimization problem, which has a better-behaved objective function. 

Remark 5: The matrix-weighted sum-MSE value l^MS'£'sum({F}, {U}, {W}, {V}) is convex with 
respect to for A; G /C if we always choose according to (|26]). Indeed, we can check that MSEt is 
convex with respect to Fg for all k,q £ IC. By construction, V^''*({F}, {U}, W^^^^) is a Hermitian 
and positive semidefinite matrix for k £ IC. Then, by definition W M S Esnmi{F} , {U} , {W} , {V}) is 
also convex with respect to F^ for A; G /C. 

IV. Algorithms 



are nonconvex 



The problems (TC), (WMSS-noPC), and (WMSS-PC) formulated in Section [HTB 
and in general are NP-hard. In this section, rather than attempting solving for their globally optimal 
solutions, we adopt an alternating minimization approach | |36} to develop iterative algorithms for finding 
their high-quality solutions. In each iteration, we alternatively fix all but one variable and determine 
the remaining variable by solving a single-variable optimization problem. The optimization problem 
in each iteration is always feasible since it has the outcome of the previous iteration as a feasible 
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point. After initialization, the algorithms are repeated until a convergent point is reached. Section IV-A 



presents the algorithm for solving (TC), which is denoted as Algorithm 1. Two algorithms for solving 



(WMSS-noFC) and (WMSS-PC) are presented in Section |rV-B] and Section |rV-q They are denoted 
as Algorithm 2 and Algorithm 3, respectively. 

A. Algorithm for Total Leakage Minimization {TC) 

There are three classes of design subproblems in Algorithm 1: i) receiver filter design, ii) relay 
processing matrix design, and iii) transmit precoder design. 

1) Receive Filter Design for (TC): We can rewrite the cost function as tr(W^ZfcWfc), where 

Zfc = Y.q=i T'k,qTl „ + Y.m=i m9k,mQl m" Sincc Wfc for A; G /C are decoupled in the cost function, 
they can be determined separately and in parallel by solving 

(r/:-Wfe):Wfc= argmin tr(X*ZfcX). 



It follows from |50| that a global optimum of {TC-Wk) is = y'tini'^k)- 



2) Relay Processing Matrix Design for {TC): We focus on determining Um for some m G by 
solving the following single-variable optimization problem {TC-\Jm) 

K K K 



mm 



k=l q=l k=l 

^k 

K K M , ^ 

fc=l q=l n=l ^ ^ 



q^k nj^m 
K K M 



i\ rs. ivi / \ 

+ EE E [yvi^yvk,n'^nnn,qni,^,^*\ 

k=l (7=1 n=l ^ ^ 



1 

qy^k n^' 
K 



S.t. 



tn X(^E'^"»,fc^m,fc + '7x,mIjVx.,„)x* J =r]v,m, (28) 
^ k=l ^ 



Where 77u,m = vT"" " En=i Ef=i tr (Un-H^.^^^^U;) - En=i 4,ntr(U„U;). Because of the 
special form of the first term in the cost function of (T£-Um), it is not straightforward to use the 
methods for the single-hop interference channel like those in | [35| to solve (T£-Um). 

We propose to transform {TC-Vm) into a more readily solvable form by introducing a new variable 
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Um = vec(Um) E C^^."^^. We define the following matrices that are independent of u^, 

K / K ^ T 



( \^ 



(29) 



■ 1 

K K M 

a2 



E E E WA.,nU„^„,,^:;,,, ) , (30) 

k=l q=l n=l ^ 



q^k nj^iri 



K 

^3,m 



K .T 

^ ^ '^m,k'^rn,k 

■ A:=l ^ 

Note that with probability one, A3 ^ is Hermitian and positive definite while Ai „ is Hermitian and 
positive semidefinite. Then, we use the following equalities, tr(ABA*C) = (vec(A))*(B-^(g)C) vec(A), 
tr(A*BA) = tr(AIA*B) = (vec(A))*(I <g) B) vec(A) and tr(AB*) = (vec(B))* vec(A) ||5T|, to 
transform both the cost function and the constraint of {TC-Um) into quadratic expressions of Um- The 
quadratically constrained quadratic program (QCQP) for designing \im is 



u 



{TC-Vim)- mill X*Ai.mX + a2^X + X*a2,m 

s.t. x*A3,mX = r/u,m- (32) 

This is a QCQP with a single equality quadratic constraint. It is nonconvex as well. 

In solving {TC-\im), we introduce a new variable Y 

([;;;{^x,m+i)x(^l,™+i) \[ follows that Y is a rank-one Hermitian positive semidefinite matrix with bottom 
right entry equal to 1. We can rewrite (TC-Um) equivalently as 






s.t. tr I I - I Y I = 7?u,„^ + 1, (33) 

= 1, (34) 

Y ^ 0, rank(Y) = 1. (35) 

While the cost function and all other constraints are convex, the rank constraint is nonconvex. This rank 
constraint is actually the main difficulty in solving {TC-\imU*^). Dropping this rank constraint, however, 
we obtain a relaxed version of (T£-UmU^), which is a convex optimization problem and also known 
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as a semidefinite relaxation (SDR) of (TC-UmU*^). Note that a complex-valued separable homogeneous 



QCQP with n constraints is guaranteed to have a global optimum with rank r < ^/n |52|. Therefore, 
having n = 1 constraints, is guaranteed to have a rank-one global optimum. The SDR of 

(TC-u.mU'^) can be solved, to any arbitrary accuracy, in a numerically reliable and efficient manner by 



readily available software packages, e.g., the convex optimization toolbox CVX |53 1. It is not guaranteed, 
however, that solving the SDR by the available software packages provides a desired rank-one global 
optimum of the SDR. Fortunately, we can construct a rank-one global optimum of the SDR from the 
resulting general-rank global optimum using the rank-reduction procedure in ||52|, which is an extension of 
the purification technique in | |54| . The key idea in each step of the procedure is to modify the eigenvalues 
of the general-rank global optima to remove the largest eigenvalue. Each step of the procedure gives us 
another global optima with the same eigenvectors but with one fewer nonzero eigenvalues. We notice that 
the last entry of the column vector obtained by the decomposition of the rank-one global optimum | |55| 
may be a complex number with modulus of 1. By multiplying the resulting column vector with the 
conjugate of its last entry, we obtain a desired column vector in the form of (u^ 1)^, which corresponds 
to another rank-one global optimum of (T/^-UmuJ^). We then use the vec^^ operator to get a globally 
optimal solution Um of (TC-Um) from the resulting u^. 

3) Transmit Precoder Design for (TC): We now focus on designing for some fc G /C by solving 
the following single-variable optimization problem 

\ ^ q=l m=l n=l ^ / 

S.t. tr(X*X) = pT,fc (36) 
tr [x* ( lil,^kl^m,k)^] = VF,k, (37) 

\ m=l J 

Where r?F,fc = pT"" " Ef=i E™=i tr (f*W;;,,W„,,F,) - =i 4,,n (u^U;;) . Note that (TC-Fk) 
is non-convex and in general is NP-hard. Since {TC-Fk) has two equality constraints, the use of the 
Lagrange multiplier method requires a more complicated 2-D search. 



Similar to Section IV-A2 we propose a method for transforming (TC-Fk) into an equivalent optimiza- 



tion problem and for solving for its global optimum. We start by defining a new variable = vec(Ffc) G 
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£^NT.kdkXi_ We also define the following matrices which are independent of 

. K M M ^ 

Bi,fc = Id.^f J] J] J]w:„,fcW;„W,,„W„,fcj, (38) 

^ q=l m=l n=l 
q^k 

. M . 

B2,k = Id.® f ^U*^,,Um,k]- (39) 

Both Bi^fc and ^2,k are Hermitian positive definite matrices. Using tr(A*BA) = (vec(A))*(I ® 
B) vec(A) |[5T|, we transform (TC-Fi^) into the following single-variable optimization problem 



{TC-fk) '■ min x*Bi ^.x 



s.t. ^*^ = P^^^, (40) 
x*B2,fcX = ?7F,fc- (41) 
Note that (TC-fk) is a complex-valued homogeneous QCQP with two equality quadratic constraints. 



Nevertheless, (TC-ik) is still nonconvex and NP-hard | |52| , |56[ . 

In solving {TC-ik), we introduce a new variable Y = xx*. Note that Y = xx* requires that Y be a 
rank-one Hermitian positive semidefinite matrix. In addition, since a*Ba = tr(Baa*) for any matrix B 
and any vector a | |57| , we obtain an equivalent optimization problem of (TC-ik) as follows 

(TC-m): min tr(Bi,fcY) 

s.t. tr(Y) = p^-^^, (42) 

tr(B2,fcY) = 77F,fc, (43) 

Y ^ 0, rank(Y) = 1. (44) 

Similar to solving (T^C-UmUj^), we adopt the SDP method for solving (T/2-f/tf^). Since (TC-fk) is a 
complex-valued separable homogeneous QCQP with n = 2 constraints, it is guaranteed that (TC-tkil) 
has a global optimum of rank r = 1 < ^/n. We can use readily available software packages, e.g., 
the convex optimization toolbox CVX |53|, to solve for a general rank global optimum of the SDR 
of (T/I-ffcf^). Next, we can always construct a rank-one global optimum of the SDR from any of its 



general-rank global optimum, e.g., by using the rank reduction procedure in |52|. The decomposition of 
the rank-one global optimum |[55| gives us the desired f^. Finally, we use the vec^^ operator to get a 
globally optimal solution of (TC-F^) from the resulting f^. 
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B. Algorithm for Sum MSE Minimization without Power Control (yVMSE-noPC) 

The design subproblems in the iterations of Algorithm 2 belong to one of the following four categories. 
1) Matrix Weight Design for {yVM.S E-noFC): Since the matrix weights V^^* for A; G /C are 



independent of each other, they can be updated in parallel based on (26). 

2) Receive Filter Design for {VV M.SS -yioFC): Recall that this approach requires the receivers use 
the linear MMSE receive filters "W^^^^ given in (pi). The receive filters can be updated in parallel. 



3) Relay Processing Matrix Design for {WM.S£-noPC): We focus on the design of Um for some 
m G 7W. By substituting (|7]) and ( [2T] ) into ( [23] ) and removing the terms independent of Um, we obtain the 
objective function of the design problem for Um- Let (WA^5<S-noPC-Um) denote the design problem 
of Um- After some manipulation and using the fact that is Hermitian, we obtain the formulation of 
(W7W5f-noPC-U„) as 

K K K 

S^h.^ ,„ E E (x^m,,^m,,X* W^^mVfcWfc.m) + 4,m E tr(X* W^,mVfc Wfc.mX) 
fc=l c/=l k=l 

K K K M 

-Y,^T{nm,kny^KrrX) + E E E tr(?^m,,^;,U;W^,„VfcWfc,mX) 

k=l k=l 5=1 n=l 



xec 



K K K M 

^tT{X*Wl^-Vk'H*^^k) + T.T.Yl tl{X*Wl^YlWk,nVn'Hn,g'H*^^,j 



k=l k=l q=l n=l 

K 



S.t. tr fx*(^?^^,fc'H;, fc + 4^^IiVx„„)x'j =r/u,m. (45) 
^ k=i ' 
Note that (WA^5£^-noPC-Um) differs from (T/I-Um) mainly due to the appearance of in the 

cost function. We introduce a new variable Um = vec(Um) and define the following matrices 

X / A' \T 

h,m ] , (46) 



Cl,m = E(E'^™''?'^™,9 + ^X,mI^x,™) ^{yVlm^kWk 
k=l ^q=l ^ 

/ K M K K X 

^2,m = vec( -J]W^,mVfc?^:;,,,+ J] J]J]W^,mVfeWMU„?^n,5^:k,J. (47) 

^ k=l n=l k=l 0=1 



Using the same manipulation as in Section IV-A2 and denoting Cs^m = A3 m, we obtain the following 
equivalent optimization problem 

(WMS£-noFC-Ufn) ■ min X*Ci,mX + C2 m^ + X*C2,m 

xec 

s.t. X*C3,mX = r/u,m- (48) 

Note that (>VA^5£^-noPC-Um) has exactly the same form as (TC-Um), thus we can apply the same 



SUBMITTED TO IEEE TRANSACTIONS ON SIGNAL PROCESSING 17 

method used for solving {TC-Um) to find a globally optimal solution of (W MS 6 -noPC-Um)- We 
then use the vec~^ operator to get a globally optimal solution of (WMSS-noPC-Vm) from the 
resulting u^. 

4) Transmit Precoder Design for (WA^iSf -noPC); We define the following matrices for the design 
of Fjfc for some k e K. 

K M M 

Dl,/c = EEE"-,fe>^UV.W,,nWn,fe, (49) 
q=l m=l n=l 

M 

D2,fe = Y.^*m,km,mn, (50) 

m=l 

M 

D3,fe = Y.^*m,kUm,k. (51) 

m=l 

After some manipulation, we obtain the following single-variable optimization problem 

{WMS£-^oPC-¥k) : min tr(X*Di feX) - tr(D^ j^X) - tr(D2 fcX*) 

s.t. tr(X*X) = pjp^^, (52) 
tr(X*D3,feX) = r?T,fe. (53) 

\ [it 1 ) = 

It follows that Y is a rank-one Hermitian positive semidefinite matrix with the bottom right entry equal 
to 1. Then, we transform (>VA1»Sf-noPC-Ffc) equivalently into the following problem 

\-(vec(D2,fc))* 1 
s.t. tr(Y) = p^'^ + 1, (54) 

(55) 

tr I I "^T,fc(ifexArT,fe<ifc QAfT.fcfifcXi j Y I _ ]^ (55) 

Y ^ 0, rank(Y) = 1. (57) 

Note that (WA^^f-noPC-ffcf^) has the same form as (T>C- u^uj^) but with one more constraint. Since 
(>V7VJ5£^-noPC-ffcf^) has n = 3 constraints (excluding the rank-one constraint), its SDR obtained 
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by relaxing the rank-one constraint is exact |52|. Thus, we can use the same steps as those in solv- 
ing (T/I-UmU^) to find the desired column vector corresponding to a rank-one global optimum 
of (VVA^5<S-noPC-ffef^). We then use the vec^^ operator to get a globally optimal solution of 
{W MS £-noFC-Fk) from the resulting f^. 

C. Algorithm for Sum MSE Minimization with Power Control (WA^5if -PC) 
In this section, we discuss briefly how to solve {yVM.S£-PC). Using the same steps as in Section 



IV-B we can develop Algorithm 3 for finding high-quality solutions of (yVMSE-VC). The details of 



Algorithm 3 ai^e provided in |38|, thus here we only compai^e and contrast the steps of Algorithm 3 and 
those of Algorithm 2. First, the matrix weight and receive filter designs for {yVM-SE-VC) are exactly 
the same as those for iyVMSS-noFC). Second, the relay processing matrix design for {WM.SS-VC) 
can be solved by the Lagrangian multiplier method with the only difference is that the multiplier must 
be nonnegative. Finally, the optimization problem for the transmit precoder design for {WM.SS-VC) is 
obtained by replacing the equality constraints in {W MS £-miPC-Y}S) by the corresponding inequality 
constraints. Fortunately, the resulting optimization problem is convex with respect to F^. In particular, 
it follows from Remark |5] that the objective function of the problem {W MS£-PC-Tk) is convex with 
respect to F^. In addition, since D3 ^ is a Hermitian and positive semidefinite matrix, then we can easily 
check that the constraints of the resulting problem are also convex with respect to F^. Thus, any available 
software package for convex optimization could be used to solve for its unique global optimum F^. 

D. Discussion 

In this section, we discuss the proposed algorithms in the following aspects: i) the convergence, ii) the 
quality of the solution, and iii) the assumption on power constraints at the relays. 

In terms of convergence, we are able to find a global optimum of the single-variables minimization 
problem in each iteration. Thus, the cost function of the original multi-variable optimization problem 
is non increasing after each iteration p6} , p7j . This guarantees that all the proposed algorithms are 
convergent. Note that the authors of ||49| adopt the alternating minimization to develop an iterative 
algorithm for solving a weighted sum-MSE minimization problem for the single -hop MIMO interference 
broadcast channel. That optimization problem has a differentiable objective function and a set of separable 
constraints in the main variables. Using the results from the general optimization |58 1, the authors of pQ] 
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are able to claim that their proposed alternating minimization algorithm converges to a stationary point 
of the corresponding weighted sum-MSE minimization problem, which is also a stationary point of the 
associated sum-rate maximization problem Nevertheless, the optimization problems in the paper, (T£), 
{yVAiSS-noPC), and {WA4SS-FC), have non-separable constraints due to the impact of transmitter 
precoders on the transmit power constraints at the relays. Thus, we are currently unable to make any 
strong claim about whether or not our proposed iterative algorithms always converge to a stationary point 
of the corresponding optimization problems. 

All the proposed algorithms are not guaranteed to reach a global optimum of the corresponding multi- 
variable optimization problem. The quality of the resulting solution depends on the initialization. One 
way to improve the performance of the proposed algorithms is to use multiple initializations, selecting 
the one with the best performance at the expenses of running time. 

The proposed algorithms in the current form are applicable only under the assumption of sum-power 
constraint at the relays. If individual power constraints at the relays are considered, we must formulate 



the corresponding optimization problems. For example, as in Section III to find high-quality solutions 
of the sum-rate maximization problem with individual power inequality constraints at the relays, we can 
formulate the following weighted sum-MSE minimization problem 

{W MS £ -PC-ind) : min WMSEsnm({F},{\J},{W},{V}) 

{F},{U},{W},{V} VL jiL j; 

s.t. pT,fc <p¥,r,A; = !,••• (58) 

Px,m<Px%,m = h--- ,M. (59) 

In some cases, we can use the same steps as in the previous sections to develop new alternative 
minimization based algorithms for solve the counterpart problems with individual power constraints 
at the relays, like {WA4SS -PC-ind). Note that the main difference between the sum power constraint 
case and the individual power constraint case is the extra constraints in the transmit precoder design. 
Specifically, with per-relay power constraints, the number of quadratic constraints of the resulting QCQP 
for the transmit precoder design is (M + 1) instead of 2 as with the sum-power constraint. The extra 
constraints make it impossible to use the same SDP method to a globally optimal solution of the transmit 
precoder design problems for individual relay power constraints when there are more than two relays, 
i.e., M > 3. Developing new methods to solve the counterpart problems to (TC) or {WAdSS-noPC) 
when M > 3 is left for future work. When power control is considered, however, we can still use the 
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same steps as in the previous sections to solve for high-quality solutions to (WA^5£^-PC-ind). Indeed, 
with power control, the single-variable optimization problems for designing the relay processing matrices 
and transmit precoders in solving {WAiSS-FC-ind) are convex |38|. Thus, we are always able to find 
their global optimum. 

V. Simulations 

This section presents Monte Carlo simulation results to investigate the average end-to-end sum-rate 
performance and to gain insights into the achieved multiplexing gains of the proposed algorithms. We 
consider only symmetric systems, which are denoted as {N^ x N'j:,d)^ + , where N-^^f^ = N^, 
-^T,fc = -^T> dk = d and Ny^^m = -^x for A; G /C, m G M^. The power values are normalized such that 
ctr j^, = cjx,m = 1, PtT = ^'x m = ^^^^ P^^"^ = for A; G /C, m G 7W. The channel reaUzations 
are flat in time and frequency. The channel coefficients are generated as i.i.d. zero-mean unit-variance 
complex Gaussian random variables. No path loss is assumed in the simulations, thus the average power 
of all cross-links on the same hop is the same. The plots are produced by averaging over 1000 random 
channel realizations. For each channel realization, the initial transceivers are chosen randomly subject to 
the power constraints at the transmitters and relays. The same initializations are used where applicable. 
Each iteration updates either one transmitter or one relay and then all the receive filters. The same order 
of relays or transmitters selected for updating is used where applicable, for example, in the comparison 
of the proposed algorithms. We use the CVX toolbox p3| to solve convex problems. 

For comparison, we consider the dedicated DF relay interference channel where one DF relay is 
dedicated to aid one and only one transmitter-receiver pair, i.e., K = M. Using equal time-sharing, the 
end-to-end achievable rate of a pair is defined as half of the minimum between the achievable rate from 
the transmitter to the associate DF relay and that from the relay to the receiver. We are interested only 
in the performance of DF relays when spatial interference alignment strategies, like those in p4| , | |49| , 
are applied on two hops. Although other interference alignment techniques, like asymmetric complex 
signaling |59], may improve the performance of DF relays, their impacts on DF relays are left for future 
work. Based on |[T2|, we derive an upper-bound on the achievable end-to-end multiplexing gain of the 
dedicated DF relay interference channel {N-^xNt^, d)^ as 0.5*min 
Individual power constraints at the relays are considered in the comparison of AF relays and DF relays. 
We also assume the transmit power at a transmitter or a relay in DF relay systems is equal to P. 



K(Ny^+N^) 

K+1 



K{N^+Ny^) 

K+1 
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1 ) Convergence: Fig. |2] illustrates the convergence behavior of the proposed algorithms. Fig. 2(a) 
provides the analysis of the sum power of post-processed leakage signals of Algorithm 1 for a random 
channel realization of the (4 x 4, 2)'^ + 4'^ system. We observe that the sum power of leakage signals 
decreases monotonically over iterations. Interestingly, the interference and the enhanced relay noise 
change their roles during the process of Algorithm 1. The interference is dominant at the beginning, 
however, it can be aligned and then cancelled quickly in a few iterations. After this point, the enhanced 
relay noise becomes dominant - its sum power is thousands times larger than the interference sum power. 
Unfortunately, given that many spatial dimensions have been devoted to deal with interference, it becomes 
challenging for Algorithm 1 to align and cancel the enhanced relay noise power. Intuitively, the enhanced 
relay noise can be thought of as a source of single-hop interference from "virtual uncoordinated relays" 
that impacts directly the receivers. Thus, we need to take into account both the interference and enhanced 
relay noise in the design of interference alignment strategies for the AF relay interference channel. Fig. 



2(b) provides the values of WMSEsum achieved by Algorithm 2 and by Algorithm 3 over iterations for a 
channel realization of the (2 x 4, l)'^-|-2^ system. We observe that WMSEsnm values for both algorithms 
are non increasing over iterations. Although the convergence speeds of the proposed algorithms are quite 
fast for these configurations, they might be slow for networks with large values of K or d. 
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(a) Total leakage power at the receivers over iterations of 
Algorithm 1 for a channel reahzation of (4 x 4, 2)^ + 4^. 




20 

Iteration index 

(b) Matrix-weighted sum-MSE values over iterations of 
Algorithm 2 and Algorithm 3 for a channel realization of 

(2 X 4, l)"* + 2". 



Fig. 2. Convergence behavior of the proposed algorithms. 



2 ) Comparison of the Proposed Algorithms: In this experiment, we simulate the average achievable 
end-to-end sum-rates for the (2 x 4, 1)^ + 2'^ system under the sum power constraint at the relays as shown 
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Fig. 3. Comparison of the average achievable end-to-end sum-rates of the proposed algorithms for the (2 x 4, 1)* + 2* system. 



in Fig. [3] We consider a sum relay power constraint. Thanks to power control, Algorithm 3 outperforms 
Algorithm 2 in this experiment. Both Algorithm 2 and Algorithm 3 outperform Algorithm 1 at low-to- 
medium SNR values because Algorithm 1 does not take into account the desired signal and noise at the 
receivers while the other do. Interestingly, at high SNR values. Algorithm 1 outperforms both Algorithm 
2 and Algorithm 3. Especially, Algorithm 1 can achieve a higher multiplexing gain than do the other. 
Zooming in on per-user achievable end-to-end rates, we find that for Algorithm 2 and Algorithm 3, some 
users have much smaller rates than do the others; they even turn off some data streams. This unfairness 
limits the maximum end-to-end multiplexing gains achievable by the two algorithms. Thus, Algorithm 1 
is more suitable than the others for investigating the maximum achievable end-to-end multiplexing gains 
of MIMO AF relay networks. 

3) Sum Power Constraints vs. Individual Power Constraints at Relays: In the previous experiments, 
we consider sum-power constraints at the relays. In this experiment, we consider the impacts of individual 
power constraints. Note that any feasible point satisfies the individual power constraints at the relays also 
satisfies the corresponding sum-power constraint. Based on the discussion in Section IV-D[ we focus 
on the case where power control is considered as it allows for the use of any number of relays. Fig. [4] 
shows the achievable end-to-end sum-rates as functions of the transmit power at a base station or a relay 
for both types of power constraints at the relays for the following three systems: i) (4 x 4, 2)^ + 4^, ii) 
(2 X 2, 1)^ + 2^, and iii) (1 x 1, 1)'^ + 2^^. We observe that Algorithm 3 for the sum-power constraint 
case slightly outperforms its counterpart algorithm for the individual power constraint case in terms of 
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maximizing average achievable end-to-end sum-rates. This gain is due to having more freedom in power 
allocation in the sum-power constraint case as relays may transmit at a higher value than the maximum 
transmit power at a relay in the individual power constraint case. This means that extra constraints 
added by the individual power constraints at the relays have little impact on the end-to-end sum-rate 
performance of the proposed algorithms. In the following experiments, we use only the counterpart 
version of Algorithm 3 that is designed for individual relay power constraints, which we refer to as 
'modified Algorithm 3'. 




5 15 25 35 45 

Transmit power at a base station or a relay [dB] 

Fig. 4. Achievable end-to-end sum-rates witli both types of power constraints at the relays for the (1 x 1, l)'' + 2^, (2 x 2, l)*+2'', 
and (4 x 4, 2)^ + 4^ systems. 

4 ) Comparison with Existing Strategies: In these experiments, we simulate several existing transceiver 
design strategies for the relay interference channel. For fair comparison, in this experiment and the 
remaining experiments, we consider the individual relay power constraints. Specifically, we simulate two 
strategies for the AF relay case. One is the AF TDMA distributed beamforming (BF), where all the 



relays help only one transmitter-receiver pair at a time (which is an extension of the design in |42| for 



multiple-antenna receivers). Another is the dedicated relay BF where each AF relay is devoted to aiding 
one and only one transmitter-receiver pair. This means that interference is ignored and we apply the 



joint source-relay design in |60|, |61| independently for the two-hop channels from the transmitters to 
their associated receivers. We also three strategies for the DF relay case that correspond to independent 
applications of single-hop strategies on two hops. The single-hop strategies include the following: i) 
selfish (SF) beamforming (i.e., each transmitter aims at maximizing the achievable rate to its associated 
receiver), ii) interference alignment strategy based on total leakage (TL) minimization p4]|, and iii) the 



SUBMITTED TO IEEE TRANSACTIONS ON SIGNAL PROCESSING 



24 



iteratively weighted MSE sum-rate (SR) maximization strategy |49|. 




15 25 35 

Transmit power at a transmitter Pj [dB] 

Fig. 5. Achievable end-to-end sum-rates for the (2 x 2, l)** + 2'^ system. 

Fig. [sjshows the results for the (2 x 2, 1)^ + 2^ system. Recall that Algorithm 3 outperforms all the other 
in all regions. It achieves an end-to-end multiplexing gain of 2 (which is equal to half of the total number 
of data streams). Note that we do not claim that this is the maximum degrees of freedom of this system. 



More complicated designs, for example those that can take advantage of symbol extensions |26|, may 
achieve higher end-to-end multiplexing gains. Unaware of interference, the dedicated relay strategies for 
both AF relay and DF relay cases achieve zero multiplexing gains. While the multiplexing gain achieved 
by the DF TL Sz TL strategy is zero, that by the DF SR & SR strategy is nonzero. The reason is 
that interference alignment is not feasible for the configuration on the two hops, interference cannot 
be completely eliminated using the TL algorithm. Although the SR algorithm is able to turn off some 
data streams, one data stream on each hop in this case, to make interference alignment feasible. Note, 
however, that it may turn off data streams of different pairs on two hops. Thus, on average the DF SR 
& SR strategy achieves an end-to-end multiplexing gain less than 1.5 (half of the number of remaining 
data streams when interference alignment is feasible). Finally, thanks to orthogonalization transmission, 
the AF TDM A distributed BF can achieve an end-to-end multiplexing gain of 0.5. 

Another experiment focuses on comparing Algorithm 3 with the the minimum-SINR maximization 



algorithm in [28|. Note that the algorithm in |28| is applicable only for single-antenna receivers. Fig. 
|6] shows the end-to-end achievable sum-rates of Algorithm 3 and the minimum-SINR maximization 
algorithm for two configurations (1 x 1, 1)'^ + 2^ and (1 x 1, 1)'^ + 2^. We observe that the end-to-end 



SUBMITTED TO IEEE TRANSACTIONS ON SIGNAL PROCESSING 



25 




a !if- - -iT^^ — n — < < < I 

I -25 -20 -15 -10 -5 5 10 15 
< Transmit power at a base station or a relay [dB] 



Fig. 6. Comparison of aciiievable end-to-end sum-rates of Algorithm 3 (siiown by tiie solid lines) and the minimum-SINR 
maximization algorithm in [28] (shown by the dashed lines) for two configurations (1 x 1, 1)'' + 2* and (1 x 1, 1)'^ + 2^. 

sum-rate performance of the minimum-SINR maximization algorithm increases at low transmit power (i.e., 
in the noise-limited regime) and saturates at high transmit power (i.e., the interference-limited regime). 
Thus, the algorithm achieves a end-to-end multiplexing gain of zero. This is reasonable since it is not 
designed specifically for interference management. Thanks to its capability of interference management. 
Algorithm 3 still achieves non-zero end-to-end multiplexing gains and provides large end-to-end sum-rate 
gains over the minimum-SINR maximization algorithm in the interference-limited regime. Our algorithm, 
however, performs worse than the the minimum-SINR maximization algorithm in the noise-limited regime 
where interference becomes a negligible issue. 

5) Maximum Achievable Multiplexing Gains: We fix = = 2 and d = 1. Fig. |7] shows the 
achievable end-to-end multiplexing gains achieved by using Algorithm 1 as a function of K for A'^x = 3 
and A^x = 5 for AF relays, DF relays, and direct transmission. We notice that with these values of 
A'x and when K is small, due to the half-duplex loss, both the AF relay and DF relay cases achieve 
lower multiplexing gains than the direct transmission. While the DF relay case cannot outperform the 
direct transmission, the AF relay case can achieve higher multiplexing gains when there are more than 6 
users. Thus, we can claim that AF relays help increase the achievable end-to-end multiplexing gains of 
interference channels. In addition, we observe that there exist upper-bounds on the achievable end-to-end 
multiplexing gains for all the simulated cases - AF relays, DF relays, and direct transmission. Theoretical 
investigation of the upper-bounds is left for future work. 
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Fig. 7. Achievable end-to-end multiplexing gains as functions of K for the (2 x 2, 1)^ + systems. 
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Fig. 8. Achievable end-to-end sum-rates of the opportunistic approach for the (2 x 2, l)'' + 2'' system with = {1, 2, 5, 10, 20}. 

6) Opportunistic Approach: The end-to-end sum-rate performance of the stationary points found by 
Algorithm 2 and Algorithm 3 depend significantly on the initializations. The opportunistic approach 
proposes to use multiple initializations and then chooses the one with the highest end-to-end sum-rates. 
Let N denote the number of random initializations. Fig. [8] shows the average end-to-end sum-rates for 
several values of N achieved by Algorithm 3 in the (2 x 2, 1)"^ + 2^ system. For this setting, at a transmit 
power of 30dB, the gain provided by the opportunistic approach over the non-opportunistic approach is 
6.4% for iV = 2, 13.2% for N = 5, 16.9% for N = 10, and 20.6% for N = 20. Note that the higher the 
value of A^, the larger the average achievable end-to-end sum-rates. Also, the additional gains obtained 
by using an extra random initialization decreases in N. Nevertheless, the benefits of this opportunistic 
approach come at the expense of longer running time. 
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VI. Conclusions and Future Work 

We developed three cooperative algorithms for joint designs of the transmitters, relays, and receivers 
of the MIMO AF relay interference chaimel. Algorithm 1 aims at miiumizing the sum power of the 
interference signals and the enhanced noise from the relays. Based on a relationship between MSB and 
mutual information, Algorithm 2 (Algorithm 3) is able to find a stationary point of the end-to-end sum-rate 
maximization problems with equality (inequality) power constraints. Simulations show that thanks to the 
consideration of the desired signal power and the noise power at the receivers. Algorithm 2 and Algorithm 
3 outperform Algorithm 1 at low-to-medium SNR. Nevertheless, they perform worse than Algorithm 1 
at high SNR due to unfairness in rate allocation among users. The multiplexing gains achievable by the 
proposed algorithms provide lower bounds on the total number of degrees of freedom in MIMO AF relay 
networks, which remains unknown. Also, the use of AF relays results in higher end-to-end multiplexing 
gains than both the use of DP relays and the direct transmission. 

A major limitation of our algorithms is that global CSI is required to implement them in their present 
form. Naturally this is challenging to achieve in a distributed system. We believe the results are still 
valuable, however, because they provide a benchmark for developing algorithms that relax the global CSI 
assumptions. Future work should focus on developing cooperative algorithms that require less overhead, 
have faster convergence speed, allow for lower implementation complexity, and account for channel 
estimation error 
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